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Abstract 

We use perturbative QCD to calculate the cross sections a L,T for the diffractive pro- 
duction of open charm (cc) from longitudinally and transversely polarised photons (of 
virtuality Q 2 ) incident at high energy (yfs) on a proton target. We study both the Q 2 
and M 2 dependence of the cross sections, where M is the invariant mass of the cc pair. 
Surprisingly, the result for a 1 ', as well as for a L , is perturbatively stable. We estimate 
higher-order corrections and find a sizeable enhancement of the cross sections. The cross 
sections depend on the square of the gluon density g (x, K 2 ), and we show that the observa- 
tion of open charm at the HERA electron-proton collider can act as a sensitive probe of the 
gluon distribution for x = (Q 2 + M 2 )/s and scale K 2 = (m 2 c + {^^(l + Q 2 /M 2 ) where the 
average quark transverse momentum squared (kj) ~ m 2 . As compared to diffractive J/tp 
production, open charm has the advantage that it is independent of the non-perturbative 
ambiguities arising from the J ftp wave function. We estimate the fraction of diffractive 
events that arise from cc production. 



1. Introduction 



The recent observation of high energy diffractive J/ ip photo- and electroproduction, 7 — > 
J/ipp, at HERA [H] has attracted a lot of interest. The principle reason is that the presence 
of the "large" scale Ml + Q 2 makes the process amenable to perturbative QCD, even for 
photoproduction (Q 2 = 0). (M^ is the mass of the J/ip vector meson.) Indeed for sufficiently 
high 7p centre-of-mass energy W the cross section for this, essentially elastic, process can be 
expressed in terms of the square of the gluon density. Thus, in principle, it seems to offer a 
particularly sensitive prove of the gluon distribution at small x. 



\|/(cc) ® T cU+p (0) ® \|/ J (cc) 




Figure 1: Schematic diagram for high energy diffractive J/ip photoproduction. The factorized 
form follows since, in the proton rest frame, the scattering of the cc system occurs over a much 
shorter timescale than the 7 — > cc fluctuation or the J / ip formation times. 



To lowest-order the 7*p — > J/ ipp amplitude can be factorized into the product of the 7 — > cc 

transition, followed by the scattering of the cc system on the proton via (colourless) two-gluon 

exchange, and finally the formation of the J/ip from the outgoing cc pair. The sequence is 

sketched in Fig. 1. The crucial observation is that at high energy the scattering on the proton 

occurs over a much shorter timescale than the 7 —>■ cc fluctuation or the J/ ip formation times. 

Moreover the two-gluon exchange amplitude can be shown to be directly proportional to the 
2 

gluon density g(x,Q ) with 

x = (M 2 + Q 2 )/W 2 , Q 2 = \(M 2 + Q 2 ). (1) 

In view of the importance of this connection, the corrections to the leading-order formula |fj, || 
have been studied [§, |J. It turns out that the major ambiguity is associated with the J /if) 
wave function. In particular there are sizeable normalization uncertainties which arise from 
allowing for the relativistic motion of the c and c quarks within the J/ip meson. Even though 
the normalization is not precisely determined, it is advocated [|J that the "shape" (or W 
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dependence) of the cross section for J/ip diffractive production at high energies can serve as a 
valuable probe of the small x behaviour of the gluon. 

From a theoretical point of view the study of diffractive open charm production has some 
advantages as compared to J/ip production. It avoids the ambiguities associated with the J ftp 
wave function and yet retains the quadratic sensitivity to the gluon distribution. Moreover, in 
contrast to J/ip, for open charm we can study the QCD behaviour as a function of M, the 
invariant mass of the cc system. In principle due to the heavy quark mass, perturbative QCD 
can predict both the cross sections a L ' T for diffractive cc production from longitudinally and 
transversely polarised photons. Indeed we find this to be the case. We compute both the Q 2 and 
M 2 dependence of a T and a L by integrating over the transverse momenta (k T ) of the produced 
c and c quarks, and over the transverse momenta (£t) of the exchanged gluons. Another feature 
is that the relevant scale at which the gluon is sampled in the diffractive production of open 
charm is 

(ml + 4) (l + ^) , (2) 

which grows with Q 2 . Besides (lowest-order) cc production we estimate higher-order QCD cor- 
rections arising from real and virtual gluon emissions. Recall that in the Drell-Yan process the 
0(as) contributions contain n 2 factors and that on resummation of these terms we obtain a 
significant enhancement of the cross section. In section 3 we will show that the 0(as) correc- 
tion to diffractive cc production also has a n 2 enhancement. 



2. The basic formulae for diffractive open charm production 

Here we study the diffractive production of a cc pair of invariant mass M from a photon of 
virtuality Q 2 at high jp cm. energy ^/s. The lowest order diagram for the process j*p — > ccp 
is shown in Fig. 2. The Bjorken x variable is given by 

x B = Q 2 /s. (3) 

As usual for a high energy diffractive process we introduce the variables 

Q 2 + M 2 

x F = , (4) 

s 

which is the longitudinal fraction of the proton energy carried by the "Pomeron" represented 
by the two-gluon exchange ladder in Fig. 2, and 

ft = Q2 — 

1 Q 2 + M 2 xjp' { ' 

The transverse momenta of the outgoing c and c quarks are denoted by ±k T , and those of the 
exchanged gluons by ±l T ; we are studying near-forward production. 
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Figure 2: Diffractive open charm production in high energy *y*p collisions, where z is the 
fraction of the energy of the photon that is carried by the charm quark. 



It is convenient to use light-cone perturbation theory (see, for example, ref. ||) and to 
express the particle four momenta in the form 

kf, = (k + ,k_,k T ) (6) 

where k± = ko ± k^. In this approach all the particles are on mass-shell, k 2 = k + k- — k 2 ^ = m 2 , 
and k + and k^ are conserved at each vertex. We choose to work in a frame in which the target 
proton is essentially at rest and where the other particles are fast with four momenta 

% = (q+, -Q 2 /q+, 0), 

k^ = (k+, m 2 T /k + , k T ), (7) 
= (£ + , £ 2 T /£ +1 tr), 

with 

2 2 , i 2 

m T = m + k T , 
where m is the mass of the charm quark. 



2.1. Factorization of the diffractive 7* — > cc cross section 

The differential cross section for the diffractive production of a cc pair of invariant mass M 



is 

72 



dM 2 dt 



t=0 A, A' 1 

where A4aa' is the amplitude for the production of a c and c (with helicities A and A') from the 
virtual photon. The 5-function arises because the mass M system is formed by c and c quarks 



with k + components zq + and (1 — z)q + . We may rewrite ([8]) in the form 



d 2 a 



dxpdt 



E 

A, A' 



M 2 + Q 2 r d 2 k T dz m 2 T 



M 2 



16vr 3 M 2 



5 z(l-z) - 



Win 



1 



M 2 / 167TS 2 



\M 



XX' 



(9) 



As mentioned in the introduction, the high-energy diffractive cc production amplitude A^aa' 
can be factorized into the light-cone wave function ip\y of the cc pair in the virtual photon and 
the (helicity-conserving) amplitude T\y for the scattering of the cc pair on the target proton. 
In analogy to ref. || we have 



Mxx'(k T ,z) 



dz' ipxx> {k' T , z) T\\i (k' T , z'\ k T , z) 



(10) 



where i/Ap occurs in the amplitude, since the cross section is the sum over the number of 
colours Nc = 3 of the charm quark. 

The factorization of A4 follows since the lifetime r 7 of the cc fluctuation of the virtual 
photon is much longer than the time of interaction with the gluons Tj. It is informative to recall 
the argument of why this is so. According to the uncertainty principle the fluctuation time 



r 7 ~ 



AE 



1 



q_ - ki_ - k 2 - 



*(1 ~ 
Q ~\~ k 2 -* 



where k\ and k 2 are the four momenta of the quarks of mass m, and 



^2 



Q 



z{l - z)Q 2 + m 2 . 



(11) 



(12) 



An estimate of the interaction time can be obtained from the typical time for the emission of 
a gluon with momentum £, from the quark ki, say. Then 



1 



jfei_ - k[_ - 



Q+ 



mtp jz — m\ I z' — l\ja 



(13) 



where a = l+/q+ and z' = z — a. In the leading logl/x approximation we have a z and 
hence 

"<?+ 



< r 7 . 



At high Q 2 the argument becomes particularly transparent. Then from ((TTJ) we have 

Q+ 1 
Tj ~ Q 2 ~ 



whereas 



gives 



Ti ~ 



m v x B 



1 

m p xi 



(14) 



(15) 



(16) 



where x^ is the Bjorken variable for the gluon-proton interaction and m p is the mass of the 
proton. We are concerned with the kinematic region < ^ where the leading log^/x^) 
approximation is appropriate, and so we have r 7 ^> T{. 
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Figure 3: Graphs which contribute to the amplitude Ty\(k' T , z'\ kx,z) for the scattering of the 
cc pair off the target proton. At high energy z = z' and the quark helicities (X and X' ) are 
conserved. 



During the short interaction time r, the exchanged gluons change only the quark (and/or 
antiquark) transverse momenta, and leave their energy fractions and helicities unchanged. Thus 
T\\i can be simply computed as the sum of the light-cone perturbation theory graphs shown in 
Fig. 3. We have § 



T xx ,{k' T ,k T ) 



Atts r d 



2N, 



c 



25{k' T -k T )-5{k' T -k T -£ T )-5{k' T -k T +£ T ) 



(17) 

where in the leading \n(l/x) approximation, z = z' and T\y becomes independent of z and z' . 
The factor arises from the propagators of the two exchanged gluons. The factor l/2Nc 
arises from the colour coupling of the gluon to the quark and the 47r from the usual coupling 
relation g 2 = Anas- The factor s arises since our amplitude T is defined so that the optical 
theorem reads "Im T = scr". The distribution f(xjp, tj) is the probability of finding a t-channel 
gluon with transverse momentum l?. That is / is the gluon density unintegrated over Gj, or, to 
be precise, It satisfies the BFKL equation, which effectively resums the leading as ln(l/x) 

contributions. We have tacitly assumed that we are considering a forward "elastic" scattering 
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amplitude with t = 0. However, the minimum value of \t\ is 

(Q 2 + M 2 V 2 2 

tmin = I ^ ™P J ~ Xjpm p . (18) 

The t min effects are expected to be small even up to xp of about 0.1 J. To relate / to the 

conventional gluon density, which satisfies GLAP evolution, we must integrate over f3p. We 
have 

xjpgix^Q 2 ) = j -^f( X]P J 2 T ) (19) 



and the inverse relation 

/(,„/•) = e T S±Ei^m. (20) 

On substituting fllTD into fllOf) we have 

M X y = i ^= J -jf Aipxx a s (f T ) f(xjp,e 2 T ), (21) 

where Aip is the combination 

Atlj(k T ,£ T ,z) = 2ij(k T ,z) - ^(k T -£ T ,z) - tlj(k T + £ Tj z). (22) 



To evaluate Atp of fl22|) we need the photon wave function [g, 0, gj. We use the convention 
of ref. || and express it in the form 



ee, \ z(l - z) 



u x (k T ,z) 7 -evyi-kr, 1 - z) 
Q ~\~ k^ 



(23) 



where ee c = \e is the charge of the charm quark, e is the polarisation vector of the photon, and 



t2 . 



Q is given by (O). 



2.2. The helicity amplitudes for diffractive cc production 

In this subsection we explicitly evaluate the helicity amplitudes which describe the 

diffractive production of a cc pair with helicities A, A' from a photon of helicity A(7). First we 
obtain the amplitudes for production from a longitudinally polarised photon and then from a 
transversely polarised photon. 

The four momentum of the photon has the form q^ = (q + , —Q 2 /q + , 0), see (0), and thus a 
longitudinally polarised photon is described by the polarisation vector 

e L = (q + /Q, Q/q + ,0). (24) 

In order to evaluate ([23|) with e = we first note that 

q^ux'jf.vx' = \ {q+uxl-vy + q-uxl+vx') = 
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so that we have 



wy_v 



Q 2 - 



(25) 



We may use this identity to evaluate the matrix element which occurs in (|23| 
u\ 'y.evy 



1 q+- Q - 

o 77 m a7-^a' H «a7+^a' 



_ 

Q+ 



2Q Jz(l-z)5 x ,-> 



(26) 



see, for example, ref. ||. Using (|23|) and fl26|) we can evaluate A?/> for longitudinally polarised 
photons. We obtain 



= -2<y A ,_v ee c Q ^(1 - 2) 



1 



1 



Q 2 + k 2 Q 2 + (k T - i T ) 2 Q 2 + (k T + i T ) 2 



(27) 



We substitute Aip L in (|21|) and carry out the angular integration. We find the amplitudes for 
the production of c and c quarks with helicities A, A' from a longitudinally polarised photon to 
be 

ee c Qz(l-z) [i ZJ- ) T\{k 2 T ) 



M L x y 



AS 



A, -A' 



/ N I 



c , 



where the integral over the transverse momenta of the exchanged gluons is 

1 1 \ 



Here we have introduced 



,Q +k 2 Q 



kt , 



Q 



(28) 



(29) 



(30) 



In order to calculate the helicity amplitudes for diffractive cc production from a transversely 
polarised incoming photon it is convenient to evaluate |7], || the matrix element U"y T • etv using 
as a basis the circular polarisation vectors of the photon 

1 



V2 



(o,o,i,±i). 



Then we obtain 

U\(-y T .e±)v X i 



:i-z) 



6 x ,-\> {(1 - 2z)X +1} e± • k T + 6 X y m\ 



(31) 



(32) 



where here quark helicities of ±~ are represented by A = ±1. In this polarisation basis we have 



A^±, = -ee c 5 Xi - X ,{(l-2z)X+l} 



-ee, 8x\> m\ 



2e± ■ e± ■ (k 

5 



T — <-T) 



_ e± ■ (kr + It) 
2 ■ k 2 T Q 2 + (k T - l T f Q 2 + (k T + £ T ) 2 



Q 2 + k 2 T Q 2 + (kr - £ T ) 2 Q 2 + {kr + £ T ) 2 



(33) 



As before we substitute Aip^ into (pi] ) and carry out the angular integration. We obtain the 
helicity amplitudes 




where the integral X\ is given by ( p9|) and 

Z 2(4 ) = / f a s (4) - ^ + ^|^) (35) 

with Q and Q\ t defined as in (0) and (|30|) respectively 

Asymptotically, when Q + k 2 -. ^> the £y integration gives a logarithmic contribution 

2 

of the form ln[(Q + fc^)//i 2 ], or rather like I/7 if the anomalous dimension 7 of the gluon 
distribution is taken into account. At realistic energies non-logarithmic contributions are ap- 
preciable and the integrations Z li2 (/c^) have to be performed explicitly. The important domain 
of integration is fi 2 < t\ < Q + k\. 



2.3. The diffract ive j TL — > cc cross sections 

To evaluate the cross sections for open charm production we need to sum |A"1| 2 over the 
quark helicities A and A'. For production from longitudinally polarised photons we have from 



E |m£v 

A, A' 



4vr 4 s 



4„2 



32 eX z\\-zYQ l I J; 



(36) 



while for transversely polarised photons we must average over the two transverse polarisation 
states A(7) = ±. From ( |34| ) we have 



1 



~ E 



A( 7 )=± A, A' 



47T 4 S 



4^2 



e 2 e 2 Ak 2 T {z 2 + {l-z) 2 }\X 2 \ 2 + m 2 \1, 



(37) 



We are now ready to calculate the cross sections from ([]). On carrying out the z integration 
in @ we find that the 6 function gives a Jacobian factor 2(1— 4m 2 -./M 2 )~2 ; and moreover implies 
that 



z(l-z) 



z 2 + (1 - z) 2 



where we have defined 



m 



r 



M 2 ' 

2m 2 T 2(3 m 2 T _ . 2/3K 2 



m 2 ' 1-/3 g 2 g 2 



K> . (38) 



s 



The variable K 2 turns out to be the scale probed by the process, since after the z integration 

Q 2 + 4 = Q 2 % + ™ 2 T = s R2 - ^ 

Collecting the factors together and inserting into @ we find that the cross section for diffractive 
open charm production from longitudinally polarised photons is 



d 2 a L 
X]P dxpdt 



An 2 e 2 c a Q 2 ^ dk 2 frr^V ^ (4Q) 



P Jo Jl-4m 2 T /M 2 \ M \ 



and from transversely polarised photons is 



2^T 



d 2 a 



dxjpdt 



7: 2 e 2 c a r\M*-m* dk ^ m \ f ( 2(3K 2 







3(1-/5) Jo ^l-Am 2 T /M 2 M 2 \ 1 \ Q 2 

(41) 

where recall that m 2 -, = m 2 + k\ and (3 = Q 2 /{Q 2 + M 2 ). The masses m and M are those of 
the c quark and the cc system respectively, and a is the electromagnetic coupling. 

2.4. Leading logarithmic approximation 

In practice the T\p, integrations over i\ of (|29|) and (|35|) have to be performed explicitly, 



> 

with the main contribution arising from the domain t\ < Q + k\. However, before we do 
this, it is informative to derive analytical expressions for the integrals assuming that the main 
contribution comes from the region i\ <C Q + k\. That is we expand the terms in brackets in 



( |29| ) and (^) in the form + c^l^ + . . . and retain only the C\ term. In this approximation 



Xi 2 are given by 



2 

T{ LA (k 2 T ) = ® ~ *f a s (Q 2 + k 2 T ) xjpgix^Q 2 + (42) 

2 

^ LA (4) = ^2^7T^« s (Q 2 + 4)^(xiP,Q 2 + 4) (43) 



where we have used ( p~9| ) to express the answer in terms of the conventional gluon distribution. 

If we use these results to evaluate \M.\ 2 and substitute them into (0) then we find, in this 
approximation, that the cross sections are given by 

d 2 a L _ 

H o ~~ 



dxpdt 

An 2 e 2 a 3 r\M 2 -m 2 dk 2 [m 2 , - 2(1 - f3)k 2 T } 2 ( - 2 rf 

-WT^tP / 1 = ■■ / 2N 3 \a s [K 2 )x w g{x P ,K 2 ^- 

3Q 4 Jo J\-Am 2 T M 2 K) 3 V 



(44) 
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k 2 (m 2 + (3k 2 T ) 2 (1-2(3K 2 /Q 2 



+ 



m 2 [m 2 T - 2{1 - f3)k 2 T } 2 ' 
4 (m|) 5 



a s (K 2 ) x IP g(x ] p,K 2 



(45) 



where the scale K 2 is given by (|38|) . 

The scale i^ 2 for diffractive processes has been emphasized by Bartels et al. and by 
Genovese et al. flOfl . Ref. |9|] concerns light quark- ant iquark production and the "hardness" 
of the scale K 2 , and the validity of perturbative QCD is ensured by considering quark jets at 
large k\. Ref. JTU] discusses open charm production, based on earlier work starting from ref. 
11|. It contains similar cross section formulae to (40) and (41), together with the logarithmic 
approximations of the form of (|44|) and (f45l). In particular the hardness of the scale K 2 = 



(m 2 + fc 2 -) /(l — j3) is emphasized, although the numerical predictions are made using the leading 
logarithmic approximation only, with the scale K 2 set equal to m 2 . Here, in section 4, we extend 
the numerical treatment to calculate explicitly the integrals X 12 over the exchanged gluon 
transverse momenta £t, and moreover we evaluate the higher order contributions (of section 
3) . We shall see that the predicted values of the cross sections for diffractive cc production are 
considerably enhanced by both of these effects. 



2.5. Connection with diffractive J/ip production 

It is instructive to see how these approximate cross section formulae given in ([P]) and 
( |45| ) compare with the result which was derived for the exclusive production of heavy quark- 
antiquark mesons [2-4]. A convenient way to make the comparison is to consider the production 
of the mesonic states with M 2 = M v <C Q 2 where 



Q 2 



Q 2 + My. 



The kinematic region for such a reaction corresponds to 

Q 2 

at large Q 2 /M v . In this kinematic region we can rewrite the cross section formulae 
(|45|) in the form 



(46) 

(47) 
and 



2^L 



M 



d 2 a 



dM 2 dt 



;i - p) x F 



d 2 a 



dxjpdt 
47r 2 e 2 c a M V Q 2 



\M 2 -m 2 



— o" \ a s{K 2 ) xjp g(xip,K 2 
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4vr 2 e 2 a M 2 Q 2 



(Q 2 + M 2 ^ 4 



2JT 



M 



dM 2 dt 



4vr 2 e 2 a M 4 



0* 



V) 
M 2 - 



«s(Q 2 ) x F g(x]p,Q 2 ] 



Tr 



(4* 



M?Am 2 + 2k 2 



T) 



4(m^) 3 



a 5 (i^ 2 ) x 1P g{x F ,K 1 



47r 2 e 2 a 



M 4 



3 (Q 2 + My) 4 



2 2 

«s(<5 ) x P g(x IP ,Q \ 



where the scale Q" is obtained by neglecting k\ in if 2 , that is Q 2 « |(Q 2 + My-), and where 
2i,,t denote the integrals over /c 2 -. For the exclusive diffractive production of a meson of mass 
My the integrals Tl,t should be replaced by the overlap integrals of the virtual photon wave 
function with light-cone wave function of the cc pair in the meson [2-4] . 

Formulae fl48Q and fl49|) have precisely the structure found for the exclusive diffractive pro- 
duction of a meson of mass My (see, for example, the J/ip diffractive production cross section 
given in eq. (2) of ref. 0]), including the result 



(49) 



olIot ~ Q 2 /M 2 . 



(50) 



However, to obtain reliable estimates of (Jl,t{J/iP) we need to evaluate the overlap integrals 
which replace T LT in ( fjgp and fl4T3|). We return now to the subject of this paper, namely the 
diffractive production of open charm. 



3. Higher order contributions to diffractive cc production 

So far we have considered only lowest-order diffractive cc production. In this section we 
evaluate higher-order QCD corrections to the cross sections. 

3.1. Real gluon emission contributions 

We wish to calculate the contributions in which the incoming photon produces a ccg three 
jet configuration. As before, the time during which the two gluon exchange interacts with the 
ccg system is much less than the lifetime of the ccg fluctuation of the photon. So we can write 
the amplitude in factorized form and consider just the coupling of the exchanged gluons to the 
ccg state. 

We need to consider only diagrams with space-like evolution. Those with time-like evolution 
do not change the cross section apart from the corrections discussed in the next subsection. 
We divide the contributions according to the transverse momenta ordering of the c, c and g. 
First we have contributions when the transverse momentum of the final gluon is much less than 
that of the quarks. The two exchanged gluons can couple to any of the three outgoing partons. 
There are eighteen such diagrams, two of which are shown in Fig. 4. These diagrams initiate 
GLAP evolution of a diffractive state starting from a gluon of transverse momentum pj*. They 
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S 1 . ,2 



} 



1 — i 3 >M 

Pt 3 Pt S 1 

Figure 4: ccg contributions with strong ordering with the gluon having much smaller transverse 
momentum than the c and c. 



give a dM 2 /M 2 behaviour for M 2 3> Q 2 ,m 2 . In contrast in the kinematic region where the 
outgoing gluon has transverse momenta greater than one of the quarks we have a dM 2 /M 4 
behaviour. Fig. 5 shows such a contribution. It corresponds to GLAP evolution of a diffractive 
state starting from an initial quark (with transverse momentum p?). Strictly speaking we 
should calculate the ccg contributions over all regions of phase space. At present the required 
formulae only exist for strong ordering of the transverse momenta. However, this should give 
a good estimate. In fact the situation is better than it first appears. For the low M 2 region 
lowest-order cc production is dominant. For higher M 2 we can use the so-called Pomeron-gluon 
splitting function, P g p, which is given by HIT 



gP{ 



IN 2 



N 2 C 



(1 + 2^) 2 (1 



(51) 



which is valid either for transverse-momentum ordering (with the outgoing gluon having the 
smallest transverse momentum) or in the BFKL limit of large \n(l/z). The latter corresponds 
to large M 2 where the ccg configuration becomes important. 

Using the formulae of ref. [0 with the coefficient functions of refs. 

d 2 a T > L 



T|, y| we find 



dxipdt 



47r 2 e 2 a a|(/i 2 



64vr 



2(3 



1 dz 



a 




xjpgixp, hp] 



where the lower limit z m is determined by the kinematical boundary 

/ Am 2 \ Q 2 + 4m 2 



01 + 



Q 2 + M 2 



(52) 



(53) 



and the cut-off Ql is chosen to be 1 GeV 2 . The leading-order coefficient functions are given by 



C J (x,r) + C L {x,r) 



x 



+ (1 - x) 2 + 4z(l - 3x)r 



3 2 2 

ix r 



In 



l+v 

1 - V 
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Figure 5: JTie ccg diagram driven by the quark distribution (jp^) 



+ v 



— 1 + 8x(l — x) — 4a; (1 — x)r 



(54) 



C L (x,r) 



4t> x(l — x) 



■^x 2 r In 



1 +v 



where r = m /Q and u, the cm. velocity of the c quarks, is given by 



1 



4m 2 x 



(55) 



(56) 



Q 2 (l-a;)' 

We take the mass factorization scale /i 2 = 4m 2 , since it was shown [15| that this is the most 



appropriate choice for the perturbative stability of (|52|) . 

From (|54D and d55|) we can gain insight into how gluon emission contributes to the evolution 
equation. For large Q 2 /m 2 we can rewrite (]54]) and ( j55|) in the form 



C T (x,r) = [x 2 + (l-x) 2 ] In 



Q 2 



•1 + 4x(l - x) 



x 2 + (1 - a;) 2 



In 



1 — x 



X 



(57) 



C L (x,r) = 4ar(l-x). 



(58) 



Since z m —> j3 in the limit of large Q 2 /m 2 , we see that the log(Q 2 /m 2 ) term in (p7|) generates 
the usual GLAP evolution equation, with the appropriate splitting function, for the diffractive 
dissociation structure function (as given in (|74|) below). 
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In principle it is straightforward to also evolve from the quark starting distribution and to 
determine its contribution at large Q 2 . However, it should be a small correction in the HERA 
regime. 



3.2. Virtual contributions 

So far we have discussed higher order contributions arising from real gluon emission. We 
should also study virtual loop corrections. For the Drell-Yan process such contributions change 
the cross section by a factor of 2 or more — the famous K factor. We must therefore investigate 
whether or not a similar enhancement occurs in the diffractive production of open charm. 
Unfortunately at present there are no complete calculations of the 0{ag) corrections for the 
diffractive process. However, it is well known that a large (usually the dominant) part of the 
K factor (say, in the Drell-Yan qq — > 7* process |16|) comes from the analytical continuation 
of the Sudakov form factor [I7H, which leads to a contribution proportional to (i7r) 2 . This 



contribution comes from the product of two imaginary parts — that is the discontinuities 
shown by the dashed lines in Fig. 6a which each give a factor of in . It is closely related to the 
double logarithmic contribution. 

At first sight it appears that there will be no it 2 enhancements to the 0{cts) corrections to 
diffractive cc production, which from a simplified viewpoint looks like the process 7* — > cc. If 
this simple view were true then to 0(as) the cross section would be a (l + as/ft), and there 
are manifestly no 7r 2 as/n terms apart of course from the usual threshold corrections. However, 
as we shall see, when we allow for the two-gluon exchange interaction we find that there are 
new diagrams which lead to a 7r 2 enhancement. 



— I — ? 1 — 



Pl I C ' P2 
I 





Figure 6: (a) Virtual and (b) real emission contributions to the Drell-Yan process q{pi)q{p2) 
7*(?)- 



Let us first recall the origin of the 0(7r 2 as/n) term in the simplest, Drell-Yan, case, which 
arise from the diagrams of Fig. 6. We first consider the emission of a single soft gluon, shown 
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in Fig. 6b, with energy uj and transverse momentum Kt- The probability to emit the gluon is 

jM 2 cIk 2 t r\ M du a s a s C F 2 M 2 



n = AC F r / — ^ = ^ In 2 ^, (59) 

J Jk, t <jJ 47T 4:71 jj, 2 

where here M is the invariant mass of the Drell-Yan pair, and where for simplicity of 

presentation we ignore, for the moment, the running of as(n T ). Such soft gluons are emitted 
independently. Summing over all possible emissions, we obtain a Poisson distribution with 
average multiplicity n. The normalization factor exp(— n) comes from the loop diagrams where 
the emitted gluon is absorbed by the other quark. The factor is [DJ 




exp 1„ » In » (60) 



where each ln(M 2 /p 2 ) contains in from the negative value of the quark virtuality p 2 < 0.Q 
Now on taking the sum of the virtual and real gluon emission terms we can obtain the 0(a>s) 
contribution to the cross section for Drell-Yan production. We have 



a = ctq 



«sCV 2 f-M 



An V // 



2 



2 



2 



where do is the lowest order cross section and /i 2 ~ \p 2 \. The factor of 2 in the last term 
arises because the real emission can occur from either the q or the q quark. Thus, since 
ln(— M 2 ) = InM 2 + in, the "real" and "virtual" «5ln 2 M 2 cancel and we have 

/ a s C F 2 
a = -o^l + ^vr +. 

= a (l + \ a s C F 7T + . 

cr exp Q a s C F 7T^ , (62) 

where the last result corresponds to the resummation of the gluon emissions. If we, as we 
should, use running q:s(k t ) in the integral of (p9|), then we obtain a ln(M 2 ) ln(ln M 2 ) form. 
Proceeding as before, and noting lnln(— M 2 ) = In InM 2 + (in/ InM 2 ), it is straightforward to 
show that the enhancement factor is such that the Drell-Yan cross section becomes 

a — (To exp(asC F 7i), (63) 

with the argument of as equal to M 2 /4. Note the absence of the factor | which occurs in the 
fixed as result fl6"2]). The larger enhancement is expected since the running of as(^ F ) m dH) 
weights the integrals towards smaller n T values. 

The same tt 2 enhancement occurs in diffractive open charm production^. It arises (in the 
Feynman gauge) from the product of the imaginary parts which correspond to the disconti- 
nuities indicated by the dashed lines in the diagrams shown in Fig. [7|. Here the exchanged 

1 There are no in terms for the annihilation process e + e~ — ► qq where the quark masses pf > 0, or for elastic 
scattering where M 2 = q 2 < and p 2 — m 2 < 0. The cut-off fi 2 in ( |59| ) is p 2 ~ \p 2 \. 
2 We are grateful to A. Kataev for discussions concerning the virtual corrections. 
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(a) 



(b) 



Figure 7: The (virtual) diagrams responsible for the n 2 enhancements to diffractive cc produc- 
tion. The exchanged gluon has virtuality £ 2 < 0. 

gluon has virtuality £ 2 < 0, whereas the quarks have p 2 > 0. Note that i 2 is the counterpart 
of M 2 for Drell-Yan production. These are the only graphs which have negative values of the 
arguments (£ 2 ' jp 2 < 0) of both the logarithms of the counterpart of (|60|) . Without the second 
t-channel gluon (indicated by £t m Fig- 0) the diagrams are analogous to ordinary deep in- 
elastic scattering for which there are no large O (asvr) enhancements. Of course there may be 
other 0(as) contributions but these are expected to be much smaller than the 7r 2 terms. Thus 
a good estimate of the cross section for diffractive open charm production is 



as in Drell-Yan production, but in this case the 7r 2 enhancement arises from quite different 
diagrams. For the diffractive process the argument of as is of the order of the largest mass 
squared of the quark-gluon state, that is m&x(p 2 ,pl), which is of the order of k\ + m 2 . The 
largest value is M 2 , where M is the invariant mass of the cc system. Here we take cts{M 2 /4). 

The 7r 2 enhancement encapsulated in (Q) , which arises from the diagrams of Fig. 7, is not 
what we would naively expect. In the space-time picture of section 2.1 we showed that the 
interaction time T{ of the cc pair with the proton is much less than the 7* — > cc fluctuation time 
r 7 . We may therefore hope that the corrections of Fig. 7 are suppressed by a factor x;/t 7 <^ 1. 
This is indeed true for the contributions coming from the real part of the diagram, but it is not 
correct for those arising from the imaginary part, that is the crucial iir terms. The imaginary 
part (the discontinuity) of the amplitude corresponds to the possibility of producing a real (not 
virtual) intermediate state which may live infinitely long. To obtain the imaginary contribution 
in5(p 2 — Trigg) we must integrate over the entire time interval from r = —00 up to r = +00 for 



a 



.cc 



o-Q C exp(a s C F Tr), 



(64) 
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the Feynman diagram in the coordinate representation. When we evaluate the discontinuity 
using the unitarity relation, 2 ImA = AA*, in the space-time picture the time in A* goes in the 
reverse direction to that in A. From the formal point of view the inverse direction of time arises 
from the opposite signs of is in the Feynman propagators (p 2 — m 2 g ± ie)~ l of the amplitudes 
A and A*. This point has been discussed in detail in the appendix of a paper by Gribov jl8| . 

Thus the large K factor of 1 + 0(7r 2 as/ir) does not contradict the space-time picture and 
the factorization properties described in section 2.1. In other words this 0(1:0.3) contribution 
may be regarded as what is left after the cancellation of the real and virtual gluon emission 
contributions which occur during the time-like evolution (or parton showering in the final state), 
which is usually not taken into account, as the leading logarithmic contributions for inclusive 
cross sections are cancelled. 

It is interesting to contrast the enhancement of diffractive open charm production with the 
situation for diffractive J/ip production. To calculate the probability for the exclusive J/ip 
process we have to convolute the cc amplitude Ai\\> of ([TI]) with the wave function of the J/ip 
meson just after the interaction. Here there is not sufficient time and energy to form a real 
intermediate quark-gluon state (analogous to the states formed in Fig. 7). Thus we do not 
expect a large K factor for diffractive J/ip production. The corrections coming from the loop 
diagrams may be treated as corrections to the J/ip wave function (including the possibility of 
ccg states in the J/ip meson). 

We now return to diffractive open charm production, the subject of this paper. We have 
considered the virtual corrections to cc production. For completeness let us also consider the 
corrections for ccg production. In this case we are less able to estimate the enhancement factor. 
In an important region of phase space the (cc)g system is produced in a colour octet-octet 
configuration, with the t-channel gluon i interacting with the s-channel gluon jet, as in Fig. f|a. 
Similar arguments lead to an enhancement of the form of (|63|) but with the colour factor Cp 
replaced by Ca 

a c ~ C9 = af 9 exp(a s C A n). (65) 

On the other hand if the colour triplet-triplet configuration c(cg) is appropriate then the en- 
hancement factor is that given in (|63|). However, the ccg contribution is only appreciable at 
the large values of M 2 , see section 4, and so the large value (and the large uncertainty) of 
the enhancement factor does not have serious impact on our main open charm predictions for 
HERA. To illustrate the ambiguity we show results for two possible enhancements of the ccg 
contribution. First, as an upper limit, we use (^). Then in a crude attempt to allow for the 
c(cg) and c(cg) configurations we take 



a CC9 = a c n cg 



1 2 

- exp(a s C A n) + - exp(a s C F ir) 



(66) 



which may be regarded lower limit. 



The exponential enhancement factor for cc production in (34) is important. It is typically 
in the range 2.7-4.0. However (Bl) is only an estimate of the true K factor. To remove 
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the scale dependence we need to know the two loop diagrams. Even at one loop level only 
the simplest, although the dominant, contribution (proportional to 7r 2 ) is taken into account. 
Nevertheless experience of the Drell-Yan process leads us to expect that the true K factor is 
well approximated by the simple expression in flBH]) . 

For diffractive J/ip production we do not have to consider the soft gluon emissions, since 
they are already included in the normalization of the J/ip wave function. The ambiguity in 
the estimates of the relativistic corrections to the wave function M is, in a sense, similar to the 
uncertainty due to the neglect of the higher loop etc. corrections to the simplified expression 
for the K factor. 

4. Numerical estimates of open charm production 

We use the QCD formalism developed in the previous two sections to estimate the rate of 
the diffractive production of open charm from both longitudinally and transversely polarised 
photons. We discuss both photo- and electro-production. 

4.1. A first glimpse of the structure of the 7* — > cc diffractive cross sections 

Before we present the detailed perturbative QCD predictions for the cross sections for diffrac- 
tive open charm production it is informative to use the leading logarithmic approximation, fl44|) 
and (|45|), to crudely anticipate some of the general features of the results. An important ad- 
vantage of the diffractive production of heavy, as opposed to light, quarks is the dominance of 
small-distance contributions. The kx integrations in (44) and (|45|) are infrared safe, protected 



by the mass m of the charm quark. Moreover the scale K 2 at which the gluon distribution is 
sampled is P, [TO] 



K 2 = ^±f (67) 

which may be considerably in excess of m 2 depending on the value of (3 and the dominant 
region of k\ sampled. Fig. 8 shows the integrands Il,t{^t) of the integrals in ([44]) and (^) for 
various values of Q 2 and (3 defined such that 

dcTjjT [l M2 ~ m2 dk 2 

xjp — = / = I L A k T) (68) 

dxjp Jo J 1 _ 4 m 2yM 2 



where the cross sections have been integrated over t assuming the form exp(— b\t\) with slopeQ 
6 = 4 GeV~ 2 [0. We comment on these plots below. 

For moderate values of (3 (0.3 < (3 < 0.8) we can obtain from the diffractive cross section 
formulae in fl4"4|) and (^) a rough estimate of the ratio o^jaT- We find 

2(m 2 T )P 3 (3{l-(3) 2 



a L : a T 



Q 4 Q 2 
2(m 2 T )Q 2 M 4 



(Q 2 + M 2 ) 3 (Q 2 + M 2 ) 



(69) 



Elastic J/ip photoproduction at HERA is observed to have a slope b — 4.0 ± 0.2 ± 0.2 GeV 2 pC 
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I L T : a) Q 2 = 12 GeV 2 , p = 0.1 



b) Q 2 = 12 GeV 2 , p = 0.4 c) Q 2 = 50 GeV 2 , p = 0.7 




xlO 



xlO 





k 2 (GeV 2 ) 



Figure 8: The k\ integrands Ii^ik 2 ), in units of GeV~ A as defined in (WB), occurring in the 
diffractive open charm cross section formulae, ($W and for xp = 5 x 10 -4 and various 

values of Q 2 and [3. MRS(A') partons fiT$ are used. 



where {rrij) is the average value of m 2 - = m + sampled by the integration in (f44|), and where 



the factor of 2 in ai comes from an approximate comparison of the integrands of ([0]) and fl45 ) . 
Eq. ( |69|) not only gives a reasonable estimate of cl/ct (up to a factor of two) but also can be 
used as a rough guide of the Q 2 and M 2 dependence of the individual cross sections — apart, 
of course, from the M 2 threshold effect coming from the rapidly expanding kinematic region of 
integration (0, \M 2 — m 2 ). From fl69|) we see that for o~l to dominate over ot we have to go 
to large Q 2 and small M 2 . For example for Q 2 = 50 GeV 2 and M 2 = 12 GeV 2 we anticipate 
from (|69D that o^joT ~ 4Q 2 m 2 /M 4 « 3 whereas if M 2 is increased to 20 GeV 2 then ai ~ o"r- 
Here we have used (k 2 -*) m m 2 , see Fig. 8(c), and hence (m\) pa 2m 2 . These rough estimates of 
< j l/°'t are in agreement with the results of the full numerical evaluation of the cross sections 
(which we will show in Fig. 10). 

To see the typical values of k\ sampled in the integrations (0) and ( f45| ) we calculate the 
average value of In k\. We show the results in Fig. 9 in the form 



(fcfi) = exp (lnfc 



as a function of (3 for our selected values of Q 2 and xp. Unless (5 is small, we see, for both ar 
and <j£, that (k 2 ^) is of the order of m 2 ; as may be expected, for example, from Figs. 8(b,c). 
Of course as we approach the kinematic bound at (3 = Q 2 /(Q 2 + 4m 2 ), which corresponds to 
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<ki> 



(GeV z ) 




0.2 0.4 0.6 0.8 




Figure 9: The average value of (kj,) for the integrals in fcftl) and for Q 2 = 12 and 50 
GeV 2 and x P = 5 x 10~ 4 . The MRS(A') gluon is used. 



M 2 = 4m 2 , the interval of integration shrinks to zero and hence (k 2 ^) — > 0. On the other hand, 
if we consider large M 2 , so that we have a large interval for the k\ integration in (fHf ), then the 
integrand II develops a different character. It is dominated by contributions near the upper 
end of the region of integration (see, for example, the dashed curve in Fig. 8(a)). It follows that 
the value of {k 2 -) is nearer the upper limit \M 2 (as shown by the behaviour of the dashed curve 
in Fig. 9 for small f3). Hence o~l samples the gluon at high scales K 2 m \M 2 . Unfortunately in 
this domain the cross section o"l is very small. 

4.2. Evaluation of the cross sections a L,T for diffractive cc production 

We calculate the cross sections for the diffractive production of open charm from the QCD 
formula of (f5]) and (|4T|) for cc production, together with the higher order contributions arising 
from ccg processes of section 3.1 and the K factor of section 3.2. The integrals 

poo Af^ 

Uk 2 T ) = / ^a 5 (4)/(^P,4)^(4,4) (70) 

JO lj> 

over the transverse momenta of the exchanged gluons are evaluated numerically, where 
i^Ki^^, k 2 -,) are the expressions given in brackets in the formulae ( |2"9"D and flSUP for T\ and 
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X 2 respectively. The contributions from the infrared region £ T < are evaluated in two 
alternative ways. We rewrite (TO) as 

d(x F g(x P ,f%,)) 

X % = a s (Q 2 )MQl) + /, d#as(4)— ^2 J - K t (£ 2 T ,k 2 T ), (71) 



where we have used (20) to express the unintegrated distribution in terms of the conventional 



gluon density. The first estimate of the infrared contribution is made using 

Ai = xjpg(x]p,Ql) Ki(Q, k%). 

We compare this result with that obtained by performing the < Bp < integration assuming 
that as(ix) g(xp,£^) = (^T/Qo) a s(Qo) 9( x P,Qo)- The two methods give essentially the same 
results. We also test the sensitivity of the predictions to variations in the choice of Qq in the 
range 0.65 to 1.5 GeV 2 . Again we find that the results are stable. 

So far we have calculated the imaginary part of the amplitude, see (|I7D . At high energy W, 
that is small xjp, our positive-signature exchange amplitude behaves as 

T oc t(xp x + (-X]p)- X ) (72) 

arising from the small xp behaviour of the gluon. Provided that A is small, A < 0.3, the 
amplitude is dominantly imaginary and the real part can be calculated as a perturbation 

ReT _ 7T 7T dlnjxjpgjxjpjp)) 

ImT ~ 2 2 d\n(l/xjp) ' 1 ' 

We calculate the correction from (|73|) and allow for a possible dependence of A on the scale £ T 
of the gluon. The contribution of the real part is included in all the predictions shown below. 

The cross sections are shown in Fig. 10 in the form 

%ip f d(3) _ 1 da* 

with i = L,T, after integration over t assuming the form exp(— b\t\) with the observed value 
of the slope parameter 6 = 4 GeV~ 2 . Eq. (|74D relates the cross sections that are shown in 
Fig. 10 to the conventional definition of F D( - 3 \ The cross sections are plotted as a function of 
the square of the invariant mass of the cc system for six different values of (x]p,Q 2 ); namely 
X]P = 5 x 10~ 4 and 5 x 10" 3 and Q 2 = 0, 12 and 50 GeV 2 . The gluon of the MRS(A') set of 
partons |L9] is used and the mass of the charm quark is taken to be m = 1.5 GeV. Throughout 
this paperQ we take the running coupling as to be such that it gives as(M^) = 0.118. The 
uncertainty Aag(Mf) = 0.005 translates into typically a ±20% error on the predicted cross 
sections. 

Some of the predictions shown in Fig. 10 are beyond the kinematic reach of the HERA 
electron-proton collider. For example if we take the maximum j*p energy for which open 
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charm can be measured to be W = 200 GeV then for xp = 5 x 10 3 we see from (HI) that 
Q 2 + M 2 < 200 GeV 2 , whereas for X]P = 5 x 10" 4 we have Q 2 + M 2 < 20 GeV 2 . 

From Fig. 10 we see that the cross section a T is generally dominant, except for low values of 
M 2 at the higher values of Q 2 . This property is clear from Fig. 11 which shows the ratio da L / da T 
as a function of M 2 . Indeed the results for the cc contributions shown in Fig. 10 quantify the 
structure that we anticipated in fl69|). We notice that the gluon emission contributions ccg (of 
section 3.1) are relatively small for M 2 < 25 GeV 2 , especially for Q 2 ^ 0. 

Fig. 12 shows the sum of all the contributions to the cross section for the diffractive pro- 
duction of open charm for three recent sets of partons MRS(A',G) jnj and GRV |^T|. The 



sensitivity to the square of the gluon density is evident. Moreover note that, as anticipated, the 
sensitivity is greater at the smaller value of xjp. The most recent HERA data for the inclusive 
proton structure function F 2 favour MRS (A') to GRV but exclude MRS(G). The predictions 
of Figs. 8-12 correspond to the choice m = 1.5 GeV for the mass of the charm quark. For 
comparison the dotted curve in Fig. 12 shows the prediction for m = 1.7 GeV for the MRS(A') 
gluon. 



individual parton sets are, of course, evolved with their respective QCD couplings. 
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Figure 10: The cross sections (xjpda L ' T /dx]p)/4ir 2 a (in GeV~ 2 ) for the diffractive production 
of open charm from longitudinally and transversely polarised photons for xp = 5 x 10~ 4 and 
5 x 10~ 3 and Q 2 = 0, 12 and 50 GeV 2 . The continuous and dashed curves correspond to o~t and 
o~l respectively. Both the cc and ccg contributions are shown. The mass of the charm quark is 
taken to be m = 1.5 GeV. The MRS(A') gluon is used. The predictions are made using 
and (51), and the "K factor" enhancements of (6^), ^6^ ) are not included. 



23 



j|qL, (cc+ccg^^T, (cc+ccg) 

x p = 5 * 10 
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Figure 11: The ratio da L /da T obtained from Fig. 10. 
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x p d0 (cc+ccg) /dx p (nb) 

x p = 5 * 10" 
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Figure 12: The cross section (xjpda/dxjp (in nb) for the diffractive production of open charm 
obtained from the QCD formulae fcjQ), and ( \5^) using MRS (A' ,G) l\Tf% (continuous, 

dashed) and GRV \2%] (dot- dashed) partons (with m = 1.5 GeV). The dotted curve corresponds 
to a charm quark mass m = 1.7 GeV for the MRS(A') gluon. The "K factor" enhancements 
of (\62) and ffih]) are included. 
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4.3. The xp dependence of the cross sections 

For diffractive production it is conventional to plot the cross section (integrated over t) 
versus ip on a log — log plot in order to see if there is a universal x^>™ dependence with n 
independent of (3 and Q 2 . The motivation for studying such a power-like dependence on xjp 
comes from the Regge inspired approach to diffractive scattering in which n is closely related 
to the (t = 0) intercept of the Pomeron. 




Figure 13: The diffractive cross section for 7*p — > Xp for large s/M x and large M\ expressed 
in triple Regge form, where Disc M ^ means the discontinuity is to be taken across the M\ cut 
of the elastic photon-Pomeron amplitude. 



Suppose we were to assume that at high energies our diffractive process is dominated by 
Pomeron exchange. Then it follows (see Fig. 13) that the cross section for large s/M x may be 
written in the form 

~ E G,«M {£) (Mlr^ (75) 

where the sum is over all Regge exchanges which contribute to photon-Pomeron elastic scatter- 
ing including the Pomeron itself. This result applies for large M|, but we can reasonably expect 
the same form with Mx replaced by M, for large values of the cc invariant mass. Recalling 
that 

■i ~ -L, m 2 = 91 

M 2 xv P 
we see that we can rewrite the cross section in the form 



M 2 = ^(1-/3) (76) 



Ff(p,Cf,t)(±) (77) 



dxjpdt Vxjp, 

where F 2 P can be regarded as the structure function of the Pomeron. The main feature of this 
Pomeron exchange model is that the power n = 2ap(t) — 1 does not depend on j3 and Q 2 and 



so can be treated flux factor [E2T 
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It is therefore informative to display our predictions for the diffractive open charm cross 
sections 

on a In F D ^ versus In xp plot to check whether a linear x P n form is obtained. The results are 
shown in Fig. 14 for the values of (3 and Q 2 that were chosen for Fig. 8. From the figure we see 
that QCD predicts that the behaviour is only approximately linear. The values of n obtained 
from the slope of the curves at selected values of xjp are shown in Table 1. We see that n 
depends on Q 2 ,/3 and also xjp. This is to be expected. For the diffractive production of heavy 
quarks we can trust the purely perturbative QCD calculation based on two-gluon exchange; we 
are not in the regime of (non-perturbative or "soft") Pomeron exchange where n = 2ap(t) — 1. 

Table 1: The exponents n of the effective x r n behaviour of F D ^ at different values of xp, Q 2 
and (3. Q 2 is given in GeV 2 . 



Xjp 


Q 2 = 12, = 0.1 


Q 2 = 12, [3 = OA 


Q 2 = 50, (3 = 0.7 


n T n L 


n T n L 


n T n L 


3 x 10" 4 
1(T 3 

3 x 1(T 3 
lO" 2 


1.40 1.31 
1.38 1.31 
1.38 1.31 
1.43 1.39 


1.48 1.52 
1.48 1.53 

1.53 1.59 

1.54 1.59 


1.53 1.56 
1.52 1.56 
1.59 1.63 
1.59 1.63 



4.4. The diffractive structure function F 2 I ^ 2; (charm) 

It is also useful to introduce the diffractive structure function integrated over xp (as well 
as t) 

F?» (f>,(?) S g- r (79) 
4ir 2 a Jxi dxjp 

which would be the charm component of the Pomeron structure function in the Pomeron 
exchange model approach. Even though we have seen that there is no basis for this model, it is 
still possible to write an evolution equation, similar to GLAP evolution, which gives the (3 and 
Q 2 dependence of F^ 2 ^ [[12], ^3f. Note that the integral over (3 of F®^ 1 gives the contribution 
of diffractive charm production to the total deep inelastic cross section and is intimately related 
to the shadowing corrections to the deep inelastic structure function. 

In Fig. 15 we show the predictions for F^ 2 ^ as a function of (3 for Q 2 = 12 and 50 GeV 2 . 
We take the limits of the integration in ( |79|) to be x\ = 0.0003 and X2 = 0.05 so as to be 



able to compare our predictions for F®^ (charm) with the experimental measurements [^4 



for the total diffractive structure function (which includes light quark production as well 
as charm). The experimental data at Q 2 = 12 GeV 2 give values of F®^ (total) in the 
range 0.15-0.2, approximately independent of (3. At Q 2 = 50 GeV 2 the measured values are 



27 




-3 -2 -3 -2 -3 -2 

10 10 10 10 10 10 

x p 



Figure 14: Plots of F D ( 3 >(xp, (3,Q 2 ), defined in ff78\), versus xjp for the values of (3 and Q 2 
used in Fig. 8. As usual the continuous (dashed) curves correspond to Ft(Fl). The values of 
the exponent n of F D ^ ~ x^ 1 are listed in Table 1. 



F®^ (total) ~ 0.3 (0.2) at (3 = 0.2 (0.65), but again compatible within the errors with no (3 
dependence. These are the values obtained by the HI collaboration El. Similar results are 



found by the ZEUS collaboration [pq| . We thus see from Fig. 15 that in this kinematic range the 
open charm contribution is predicted to be approximately 25-30% of the measured diffractive 
structure function. The results shown in Fig. 15 are obtained using MRS (A') partons and the 
mass of the charm quark m = 1.5 GeV. If a mass m = 1.7 GeV were to be used then the 
height of the peaks shown in Fig. 15 would be reduced by 20%. If on the other hand GRV 
or MRS(G) partons are used then the peak value of F®^ (charm) at Q 2 = 50 GeV 2 is found 
to be 0.17 or 0.24 respectively, as compared to 0.074 for MRS(A'). When compared with the 
experimental values of F®^ (total), these high values of F®^ 1 (charm) clearly disfavour the 
small x gluon distributions of these sets of partons, even allowing for the uncertainty in the K 
factor enhancement (|64T) or due to the choice of m. 



Fig. 15 displays clearly the characteristic (3 dependence of diffractive open charm production. 
First we see the kinematic bound for cc production. It corresponds to an upper limit on /3, 

Q 2 

^ max = Q 2 + M 2 h (80) 

where the threshold mass is given by M 2 h = 4m 2 . The peak seen just below this value of (3 
arises from cc production just above threshold with a fall off as M 2 increases (i.e. j3 decreases). 
The rise which occurs for smaller (3 corresponds to ccg production. Recall that for the ccg 
component there is a large enhancement arising from virtual corrections to this process. The 
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Figure 15: The continuous and broken curves are the predictions for the contribution of 
diffractive open charm production to the structure function F^ 2 \f3,Q 2 ) of (I73j as a function 
of P for Q 2 = 12 and 50 GeV 2 respectively. The MRS(A') gluon is used, xp is integrated over 
the interval (t).0003, 0.05,). The curves branch at low values of (3 according to whether we use 
the enhancement factor ME) or ffih}) arising from the virtual corrections to ccg production. 



size of the enhancement is not well known. This uncertainty in the ccg contribution is repre- 
sented in Fig. 15 by two curves which correspond to two different choices, (|65|) and ( |66|) , of 
the ccg enhancement K factor. Indeed at the lowest value of (3 shown (for each Q 2 ), after the 
cc contribution is subtracted^, we find that the residual ccg contribution has approximately a 
factor of two uncertainty. 



5 The cc contribution vanishes linearly with (3 as (3 — > 0. 
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5. Conclusions 

We have presented the QCD predictions for the cross sections (<Jl,t) fc> r the diffractive 
production of open charm from longitudinally and transversely polarised photons of virtuality 
Q 2 . At lowest order, the diffractive processes j*P ~^ °cp are driven by two-gluon exchange 
between the cc pair and the proton. The perturbative predictions are protected by the mass 
m of the charm quark and are infrared safe. In fact the diffractive production of open charm 
depends on the square of the gluon density xjpg(x]p, K 2 ), with xp = (Q 2 + M 2 ) / s at a scale 

K 2 = (<4> + m 2 ) (l + (81) 

where typically (kj) ~ m 2 . M is the invariant mass of the cc system. 

We have aimed to make our diffractive cross section predictions as realistic as possible for 
the experiments at HERA. We have studied the Q 2 and M 2 dependence of both a L and <j T . 
We found that a L only exceeds a T in the kinematic region of low M 2 and high Q 2 . There 
are several new features incorporated in our analysis. First, we integrate explicitly over the 
transverse momenta ±£t of the exchanged gluons. This proves to be important, since we 
find that the inclusion of the l\ effects increases the cross sections by about a factor of 2 as 
compared to using just the leading logarithmic approximation. Second, we have estimated the 
higher-order contributions to diffractive open charm production. These were divided into the 
computation of real emission ccg contributions, and the estimation of the enhancement due 
to virtual corrections to cc production. The real gluon emission contributions were found to 
be small for low M 2 , M 2 < 25 GeV 2 , and moreover decrease with increasing Q 2 . The virtual 
corrections on the other hand, are surprisingly important. They arise from the diagrams of 
Fig. 7, and have some similarities to the ir 2 enhancement of the 0(a s ) corrections that is well- 
known in Drell-Yan production. However, their occurrence in diffractive 7* — > cc production 
is a novel and theoretically interesting effect. To a good approximation they can be resummed 
and represented by a K factor enhancement in the form of an exponential with argument 
0(as7f)- Just as in Drell-Yan production, we estimate an enhancement of the lowest-order 
result by a factor of about 3. Due to the larger colour coupling of the gluon, the enhancement 
is expected to be greater for the ccg contribution. However, in this case the K factor is not 
well-known. Thus in the large M 2 (and small Q 2 ) region, where the ccg contribution is not 
negligible, our cross section predictions have a larger uncertainty (as represented, for example, 
by the branching of the curves shown in Fig. 15 at the lower value of (3). 

We have used the above perturbative QCD formalism to calculate the cross sections <tl,t 
for diffractive open charm production. We have presented representative numerical results 
to illustrate the Q 2 , M 2 and xp dependence of the cross sections, which are relevant to the 
experiments at HERA. In particular we show the sensitivity to the choice of gluon distribution 
at small x. 

We may compare diffractive open charm production with diffractive J/ip production at 
HERA. Both processes are special in that they depend quadratically on the gluon distribution 
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at small x. The main uncertainty in the calculation of the cross section for diffractive J/i/j 
production was found to be associated with the Fermi motion of the c and c within the J/i/j 
and the choice of the mass m of the charm quark [|], [5]. The uncertainty is greater in the 
predicted size of the cross section than in the energy (or equivalently the x) dependence of the 
cross section. That is the shape, rather than the normalization, of the observed cross section 
is a better discriminator between the various gluon distributions. The shape of the diffractive 
J/ip photoproduction data collected at HERA favours the MRS (A') gluon relative to that of 
GRV, and rules out the MRS(G) gluon §. 

Diffractive open charm production has the advantage that it is independent of the uncer- 
tainties due to the J/i/j wave function. On the other hand the higher order corrections to open 
charm have (novel) it 2 enhancements which lead to a large K factor which, at present, can only 
be estimated. The K factor uncertainty is not expected to be present in the J/i/j prediction; 
it is automatically subsummed by using the observed leptonic width to fix the J/i/j coupling. 
High energy data for the two processes will therefore act as complementary and independent 
probes of the gluon at small x. Their quadratic sensitivity to the gluon means that valuable 
information can already be obtained despite the above uncertainties. 

When we compared our perturbative QCD predictions for diffractive open charm production 
with the inclusive data for diffractive production, at a given f3 and Q 2 , we estimate that about 
25-30% of diffractive events arise from cc production (if the MRS (A') gluon is used). The j3 
(and Q 2 ) dependence of diffractive cc production has special features that are well-illustrated 
in the sample results shown in Fig. 15. Low mass cc production leads to a characteristic peak 
in the region (3 < Q 2 /{Q 2 + 4m 2 ), while ccg production only becomes important at much lower 
(3. To identify the dramatic threshold behaviour it will be necessary to reconstruct the mass of 
the cc system. 

In summary we have used perturbative QCD to predict diffractive open charm production 
at HERA as a function of Q 2 , M 2 and xjp. The main unknowns are (i) the gluon distribution 
at small xjp, (ii) the mass m of the charm quark and (iii) the accuracy of the estimate of the 
large K factor enhancement. The quadratic sensitivity to (i) means that information on the 
gluon can still be obtained despite the ambiguities arising from (ii) and (iii). If the gluon is de- 
termined at small x by independent means, then the characteristic Q 2 , M 2 and xjp dependence 
of diffractive cc production will offer a valuable probe of the validity of perturbative QCD at 
xp < 10~ 3 and at scales of the order of a few times m 2 . 
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